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REFLECTION AND TRANSMISSION OF ACOUSTIC WAVES 
FROM A MOVING LAYER 


By George G. Steinmetz and Jag J. Singh 
Langley Research Center 

SUMMARY 

The refraction of acoustic waves by a moving medium layer is theoretically treated 
and the expressions for reflection and transmission coefficients are determined. The 
moving medium layer velocity is assumed to have a space dependence in one direction. A 
partitioning of the moving medium layer into constant-velocity sublayers is introduced and 
the number of sublayers is allowed to increase until the reflection and transmission coef- 
ficients converge to their respective values. 

The introduction of the input impedance concept and the use of several dummy vari- 
ables molds the equations into a form very suitable for computation. Numerical results 
have been obtained for several sublayer approximations of Poiseuille’s flow as functions 
of the moving layer velocity for several angles of incidence of the acoustic wave. The 
refined approximations are shown to converge. The degenerate case of a single constant - 
velocity layer is also treated, both theoretically and by a numerical analysis. 

INTRODUCTION 


The investigation of reflection and transmission coefficients of an acoustic wave 
from a moving layer has been treated by a number of authors. (See refs. 1 to 5.) These 
calculations indicate that the reflection and transmission coefficients are strong functions 
of the velocity of the medium and the angle of incidence of the acoustic waves. However, 
most of these treatments have been confined to medium velocities which are constant; 
that is, the medium flow is both time and space independent: 


dY=0 

dt 


ft-0 

dy 


^ = Constant 
dx 


where V is velocity, t is time, and x and y are coordinates. In many practical 
cases, the medium flow exhibits a space dependence. For instance, the liquid flow through 


cylindrical pipes has a strong velocity gradient along one space axis. This well-known 
case of Poiseuille's flow can be described by a second-order polynomial. A problem of 
this type may be approximately treated by partitioning the moving layer into a number of 
sublayers and considering the velocity in each sublayer to be constant. 

It is the purpose of this report to define suitable computational techniques for study- 
ing the propagation of an acoustic wave through a medium moving with space -dependent 
velocity. The approach is, as indicated previously, to divide the moving medium into 
constant -velocity sublayers and use the existing theoretical solutions to propagate the 
wave across each of the sublayers. This procedure involves matching boundary condi- 
tions deduced from the requirements that the pressure must be continuous and the instan- 
taneous component of the resultant velocity (stream and acoustic) locally normal to the 
rippled interface between the adjacent layers must be the same on both sides (ref. 6). 

The mathematical development is for a general n-sublayer case. Numerical results are 
given for a small number of sublayers (5) and for a large number of sublayers (25 to 45). 
An optimum number of sublayers based on convergence in the magnitudes of the reflection 
and transmission coefficients, as the partitioning is refined, are sought. A preliminary 
account of this work has been presented in reference 7. 

SYMBOLS 

A n amplitude of transmitted wave 

B n amplitude of reflected wave 

C speed of sound 

d depth of entire medium 

d n depth of nth sublayer 

i=fl 

kn = |^n| = (&> - V n k n>x )^C n 
M number of layers 

n number of sublayers, n = 1, 2, 3, . . ., M 

p acoustic pressure 
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R reflection coefficient of moving layer, nondimensional 



T transmission coefficient through moving layer, nondimensional 

t time 


*n = 


k n,y 

k M,y 


U x 

u y 

v 


x,y 

z 

z in 


cos 9 


M 


V„ 

1 - —2. sin 9 
\ c n 


M, 


- sin 2 0 



x-component of acoustic velocity 
y -component of acoustic velocity 
medium flow velocity 
Cartesian coordinates 
characteristic impedance of medium 
input impedance of medium 



Z* n = zin 
pu>d 


a 

9 

P 

cr 

U) 


modulus of y-component of propagation vector 

angle at which incident acoustic wave strikes medium (or leaves medium) 
density of fluid 

modulus of x-component of propagation vector 
acoustic wave frequency 
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‘ PROBLEM APPROACH 


Figure 1 illustrates the three -layer basic problem. Layers 1 and 3 are stationary, 
and layer 2 has space -dependent flow. Layer 2 is partitioned into M - 2 sublayers, 
each with a constant velocity. (Details of the partitioning technique are presented in a 
subsequent section.) For individual sublayers, the problem can be treated in the manner 
of the constant -velocity cases (ref. 2). Basic equations for the acoustic propagation for 
constant -velocity layers are (ref. 1) as follows: 


ap„ 9 P„ o 
^♦V nl J + c2pl 


'9U-* 


au. 


y,n 


ax 


9 y 


8U x,n 

, v 

9 U x ,n 

. 1 

a P n 

at 

+ v n 

ax 

Pn 

ax 

9U y,n 

+ V n 

9U y,n 

+ ± 

9 Pn 

at 

T n 

ax 

P n 

9 y 


(1) 

( 2 ) 

(3) 


where subscript n denotes a particular sublayer (n = M and 1 represent the top and 
bottom stationary layers, respectively). 

The general solution in a layer n = 1, 2, . . ., M is as follows: 


p n = |^n exp(-io! n y) + B n exp(ia n yjjexp|i(wt - cr n x) 
Uy, n = [V»n(* " CT n v n)|( 9 Pn/ 9 y) 


(4) 

(5a) 


or 


U- 


a 


y,n = 


n 


P n ( w - CT n v n)L 


A n exp(-ia n y) - B n exp(ia n y) exp i(wt - cr n x) 


(5b) 


where A n and B n represent the amplitudes of the transmitted and the reflected waves, 
respectively. Note that for n = 1, Bj = 0 since there is no reflected wave in that layer. 
The propagation vectors a n and a n are defined as follows: 


“n ~ k n? y - 
a n = k n,x = 



| cos 0 n 
sin & n 


(6a) 

(6b) 
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where 


v n 


= ( k n,x + k n,y) - (u - V n k njX ^C n 


(6c) 


The boundary conditions that must be met at each sublayer interface are continuity 
of the acoustic pressure 


P n+ 1 = P n < 7a > 

and the instantaneous component of the resultant velocity locally normal to the applied 
interface 


u y,n+l j(^ ~ CT n+l^n+l) “ u y,n j(^° " CT n^n) ( 7b ) 

Let each layer impedance be denoted by Z n , where 

/ tt- \ 2 2o2 

_ Pnl^ - v n^n) _ P n M Sn (R , 

n “n ‘ «n W 

with 

*2 = (' ' § 8tn 9 “) 

After applying the foregoing boundary conditions, the input impedance Z^ n is given by 
the following expression (ref. 5): 

^n = (^n-1 “ tan ot n dj^jl(z n - iZ n _i tan a n d n )z n (9) 

■in 

where d n is the nth sublayer depth and Zj = Zj. 

The reflection (at the surface of the moving medium) and transmission (through the 
moving medium) coefficients are then given by 

R - fen - Z n+ l)/(z!. n + Zn+l) (10) 
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and 


T - 1 1 [ z i + z f)/( z J+ i - <“> 

It should be emphasized that velocity dispersion and other frequency -dependent effects 
have been ignored in these calculations. 

ANALYSIS AND RESULTS 


Let the acoustic source be located in the Mth layer and let d be the width of the 
moving medium. The value of kj^d, even though M is a variable (i.e., a function of the 
number of partitioned layers chosen to approximate the moving medium), is maintained 
constant for all values of M. (See ref. 5.) For notational convenience, a new variable 
t n is defined as 


: _ kn >y _ 1 

n k M,y cos e M 



(n = 1, 2, 3, . . ., M) (12) 


The variable t n is a complex number for selected ranges of V n /C n and it is either 
real or imaginary with the passage between the two regions yielding t n = 0; t n is imag- 
inary if 

1 - sin d M < Vn < sin % + 1 
sin 0 m C n sin 0^ 

and t n is zero for 

= (1 ± sin 0 M)/ Sin e M 
Otherwise, t n is real. 

A new impedance Z n is now defined in terms of t n as follows: 

_ Z n (»n /f)4 

n w 2 pd t n k M < ^ cos % 
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The new input impedance is given by 


-in COS [ tnA ( dn / d )] ' i ( 1 / A ) s ?(P/P n )( 1 / t n) sin ^n A ( d n/ d ) J 


Z “* _ 

n - 


S n cos t n A(d n /d)^ - lZ^A^p/p^sin t^c^/d) 


(13) 


where A = k^jd cos 0^. 

If t n = 0, then limits must be taken to find the value of Z^ . The appropriate 
limits yield 



_ 

- ^n-l 


- »(P n /P) (■%/<#§) 


(14) 


Now, the resulting reflection and transmission coefficients are given by 


and 


Z M-1 ~ ( p m/ P )^ 1//A ^ 

*m-i + (p m / p ) (1/a) 


j=M-l 

T= |T — 

j=l o2 


•(Pj/P) + 


77in 


"i 2 ) 


TTin 


exp 


s j+i( p j+i / p ) + At j+i z j 


Mf 


(15) 


(16) 


Assume identical media in all layers, that is, Pj = p 2 = . . . = p M and 

Ci = C 2 = . . . = C M . The case of a single layer with constant velocity can be 

treated as a special case of the given formula. 

The reflection coefficient for the special case can be shown to vanish at points 
of V/C for fixed incident angles 9 by 



sin 0 M 


r 

2 

1/2" 

■ 

(k^d) +sin20 M 

> 

J 


a = 1, 2, 3, . . .) (17) 


If l = 0, then t n = 0 and the limit process given by equation (14) must be used. Fig- 
ures 2 and 3 illustrate the calculated values of the magnitudes of the reflection and trans- 
mission coefficients as functions of V/C for two values of 9. 


7 



Let the moving medium velocity now have a one -dimensional space dependence. 
Chosen for this example is Poiseuille's flow in which the velocity is described by the fol- 
lowing parabolic curve: 

® = y(i-y) (osysi; -2^^) (18) 

c r/max x ' 

Let the moving medium be partitioned into three equal (dj = 6.2 = d 0 ^ sublayers (M = 5) and 
the velocity in each sublayer be treated as a constant determined by the value at the mid- 
point of the layer. The geometry of this situation is shown in figure 4. Calculated val- 
ues of the magnitudes of the reflection and transmission coefficients for the M = 5 case 
are shown in figures 5 and 6 , respectively. 

The three -equal -sublayer case (M = 5) is at best a crude approximation of the para- 
bolic flow. Clearly, increasing the number of sublayers gives a better approximation. It 
appears more appropriate, however, to maintain some limiting difference between flow 
velocity in adjacent layers instead of a constant partition width, as has been done so far. 
With this new technique, a variable width partitioning was used, based on the criterion that 
the flow velocity difference between adjacent layers will not exceed a certain fraction of 
the peak flow velocity. 

Two examples of this calculational refinement are presented. In the first example, 
the adjacent -layer velocity difference was restricted to a value =10 percent of the maxi- 
mum velocity. The partitioning scheme is illustrated in figure 7. Note that the flow 
symmetry about the midstream point is preserved. The magnitudes of the reflection and 
transmission coefficients calculated for M = 25 are shown in figures 8 and 9, respec- 
tively. In the second example, the partitioning was refined even further by requiring that 
the adjacent -layer flow rates differ by =5 percent of the maximum flow rate. The calcu- 
lated magnitudes of the reflection and transmission coefficients for the M = 45 case are 
illustrated in figures 10 and 11, respectively. A comparison of the results of these two 
examples reveals that they are almost the same, which indicates a convergence in the 
respective values of the reflection and transmission coefficients. It would thus appear 
that a criterion of AV H 0. lV max and M = 25 gives a sufficiently accurate represen- 
tation of the actual flow. 


DISCUSSION OF RESULTS 

For the degenerate case of a single velocity layer, the magnitudes of the reflection 
and transmission coefficients given by equations (15) and (16) are the same as those 
reported in the literature (ref. 2). The graphical representations of these equations given 
in figures 2 and 3 agree with the zeros of |r| as defined by equation (17). 
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The pattern of the magnitudes of the reflection and transmission coefficients exhib- 
ited in figures 5, 6, 8, 9, 10, and 11 for the cases of a partitioned y -dependent flow shows 
a generally similar, but distinctly different, signature from that of the case of a single 
constant -velocity layer. The position of zeros of |r| has been altered somewhat as M 
is increased. For example, the zero value of [r| has shifted outward to higher numer- 
ical values of [v/C | for higher M values, with convergence in the highest M values. 
Equation (17) is not valid for any case but M = 3 and the method of solving for zeros 
of |r| becomes very cumbersome for more than a single moving layer. A comparison 
of the results for refined partitioning with the results for the crude three -sublayer case 
indicates very little change in contrast to the degree of difference between the results for 
the single layer and multilayer cases. The difference between the two highly refined par- 
titioned cases (M * 25 and M = 45) is hardly detectable, which indicates that the reflec- 
tion and transmission coefficients have converged to their respective values at M = 25. 

The computer solution requirements for the final equation are minimal as a result 
of the form in which these equations were cast. The results presented for each case in 
this publication represent under 45 seconds of central processing time on the Control 
Data 6600 computer system. Note that this processing includes three different angles of 
incidence and a sweep of V/C between ±2.0 with AV/C of 0.01. The processing of the 
M = 45 case required a threefold increase in time over that of the M = 5 case. 

CONCLUDING REMARKS 

The reflection and transmission of acoustic waves by a medium moving with 
Poiseuille's flow have been treated theoretically. The moving medium was subdivided 
into a number of sublayers, with each sublayer having a constant velocity. The number 
of sublayers was allowed to increase until a saturation was achieved in the values of 
reflection and transmission coefficients. 

With the introduction of the input impedance concept and the use of several dummy 
variables, the equations have been altered into a form very suitable for computer analysis. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., February 8, 1972. 
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v/c 

Figure 3.- Absolute values of the transmission coefficient as a function of the velocity of the moving medium 
for two angles of incidence. M = 3; kj^d = 2.0; p = 1.000; = 10.0°; Q<^ - 30.0°. 
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Figure 7.- Uniform velocity multilayer approximation of Poiseuille's flow (M = 25). Velocity difference between 
any two adjacent layers equals 10 percent of the maximum fluid flow rate. 
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Figure 9. - Absolute values of the transmission coefficient as a function of the velocity of the moving medium 
for three angles of incidence. M = 25; kj^d = 2.0; p = 1.000; = 10.0°; 0 2 = 30.0°; 9 g = 45.0°. 
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Figure 10.- Absolute values of the reflection coefficient as a function of the velocity of the moving medium 
for three angles of incidence. M = 45; k^d = 2.0; p = 1.000; = 10.0°; 0 2 = 30.0°; 0 3 = 45.0°. 
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Figure 11.- Absolute values of the transmission coefficient as a function of the velocity of the moving medium 
for three angles of incidence. M = 45; = 2.0; p = 1.000; 0 1 = 10.0°; d 2 = 30.0°; = 45.0°. 
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